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Abstract 

We discuss a possibility that the Neutron Electric Dipole Moment (NEDM) can be calculated 
in lattice QCD simulations in the presence of the CP violating 6 term. In this paper we measure 



the energy difference between spin-up and spin-down states of the neutron in the presence of an 



^ ■ uniform and static external electric field. We first test this method in quenched QCD with the RG 

improved gauge action on a 16 3 x 32 lattice at a -1 ~ 2 GeV, employing two different lattice fermion 
formulations, the domain-wall fermion and the clover fermion for quarks, at relatively heavy quark 
mass (mps/my ~ 0.85). We obtain non-zero values of NEDM from calculations with both fermion 
formulations. We next consider some systematic uncertainties of our method for NEDM, using 
24 3 x 32 lattice at the same lattice spacing only with the clover fermion. We finally investigate 
the quark mass dependence of NEDM and observe a non-vanishing behavior of NEDM toward 
the chiral limit. We interpret this behavior as a manifestation of the pathology in the quenched 
approximation . 

PACS numbers: 11.30.Er, 11.30.Rd, 12.39.Fe, 12.38. Gc 



1 



I. INTRODUCTION 



Discrete symmetries, such as parity (P), charge conjugation (C) and time-reversal (T), 
have played important roles to establish the structure of the standard model. One of the most 
famous examples is CP violation which led to three generations of quarks and leptons[l|. 

In the strong interaction, the most strict constraint on violation of P and T symmetries 
comes from the measurement of the electric dipole moment (EDM) for neutron (NEDM) 
and proton (PEDM). The current upper bound is given by 

\d N \ < 6.3 x 1(T 13 e ■ fm (90%C.L.) (1) 

for neutron from a Larmor frequency measurement with ultra-cold neutron (UCN)^J, and 

\d N \ < 5.4 x 10" 11 e-fm (2) 



for proton|3j, which is estimated indirectly from EDM of mercury atom 199 H g given by 
4tom( 199 H g ) < 2.1 x 10- 15 e-fm (95% C.L.)[4|. 

On the other hand, QCD allows a gauge invariant, renormalizable, and CP odd 9 term, 

^^T~? [ d A xG llu (x)G IMU (x), G^uix) = -e llua pG a p{x) (3) 

OZ7T z J 4 

in Euclidean space-time with G M „ which is the field strength of gluon. Some model 
estimations [3, S] yield 

\d N \ ~ 9 x 0(l(r 2 ~ l(r 3 )e-fm, (4) 

which leads to a bound 9 < O(10~ 10 ). Hence 9 must be very small or even may vanish in 
QCD. 

Smallness of 9 in the QCD sector, however, is not protected in the presence of the elec- 
troweak sector of the standard model, where the quark mass matrix, arising from Yukawa 
couplings to the Higgs field, may be written as 

ffWMijtfix) + #(x)M^f (x), (5) 

where ip L and ip R represent left and right handed quark fields with flavor indices Diag- 
onalizing the mass matrix and making it real, the parameter 9 becomes 

# = #QCD + argdetM, (6) 
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where #qcd is the original 9 parameter in QCD. Therefore, #qcd and argdet M contributions 
have to cancel out to the precise degree that the stringent experimental upper bound on 
NEDM is satisfied. In either of the two cases, it seems necessary to explain why Nature 
chooses such a small value for 9; this is the "strong CP problem". One of the most attractive 
explanations proposed so far is the Peccei-Quinn mechanism jf]]. Unfortunately, the axion, a 
new particle predicted by this mechanism, has not been experimentally observed so far. 
Present theoretical estimations of NEDM vary in mag nitude among different models such 
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as current algebra^], chiral perturbation theory 

also [3]). While these crude estimations of \djsr\/9 already convince the smallness of 9, a 



ll|, and QCD sum rule|12l. 1 131] (see 



theoretically reliable and accurate estimation for NEDM will be required to determine the 
value of 9, if a non-zero value of NEDM is observed in future experiments. Lattice QCD 
calculations provide a first-principle method for this task. Indeed more than fifteen years ago, 
the first attempt was made to estimate NEDM in a quenched lattice QCD simulation 15| . 
Reliable signal of NEDM could not be obtained at this time 16| . Since then, no lattice 
calculation of NEDM have been attempted until recently. In the last year, new approach has 



been presented for this problem. Ref. 13, ll8|, ll9[ proposed a formulation to extract the CP 



odd electromagnetic form factor of nucleon from certain lattice correlation functions. NEDM 
can be extracted from this form factor in the zero momentum transfer limit. Applying this 
formulation in a quenched calculation with domain-wall quarks, a non-zero value for the CP- 
odd form factor of nucleon was obtained at one value of non-zero momentum. Based on this 
formulation, the same form factor has been calculated on gauge co nfig urations generated by 
Nf = 2 dynamical domain- wall QCD at several non-zero momenta 20j]. The value of NEDM 
after the zero momentum extrapolation, however, is consistent with zero within the large 
statistical error in this calculation. 

The results mentioned above suggest that, while it is possible to obtain signals for the 
CP odd form factor at fixed and small value of momentum, it is numerically difficult to 
carry out a statistically controlled extrapolation of the form factor to the zero momentum 
limit to extract the value of EDM. Therefore, in this paper, we investigate another method 
to calculate the value of EDM directly without momentum extrapolation. In this method, 
introducing a constant uniform electric field E, we measure the energy difference between 



spin-up and spin-down components of the nucleon in the presence of the 9 term 



15|. If the 



electric field is small enough, the leading contribution to the energy difference is given by 
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— * — * — * —> 

d^S ■ E with neutron spin S and electric field E. Therefore EDM can be directly extracted 
without momentum extrapolation. The most difficult part of this calculation is to reweight 
the nucleon propagator on a given gauge configuration with the factor e l6 ®, where Q is the 
topological charge of the configuration. We may control this reweighting by taking a small 
value of 9. Another difficulty is that our electric field breaks periodicity in the time direction, 
generating large field at the time boundary. We should investigate influences of the large 
electric field at the boundary to EDM signals. 

We check the ability of this method in the quenched approximation at a heavy quark 
mass. We employ two fermion formulations, domain-wall fermion having chiral symmetry 
and clover fermion with explicitly broken chiral symmetry, in order to investigate possi- 
ble dependence of EDM signals on the aspect of chiral symmetry of fermion formulations. 
Our study have revealed that the quality of EDM signals is not very sensitive to fermion 
formulations. Therefore we have employed the clover fermion, which requires much less com- 
putational cost than the domain-wall fermion, to study various systematics of EDM such 
as the volume dependence, the boundary effect and the quark mass dependence within the 
quenched approximation. 

This paper is organized as following. In sec. [Ill we explain the definition of EDM and 
our method to extract EDM from nucleon propagators. Simulation details of our lattice 
calculation are summarized in sec. IIHI In sec. IIVI we show numerical results with both 
domain- wall and clover fermion at heavy quark mass on a 16 3 x32 lattice. We then investigate 
the finite size effect and the boundary effect on a 24 3 x 32 lattice with the clover fermion. 
In [V] we systematically study the quark mass dependence of EDM using the larger lattice 
with the clover fermion. A summary and discussion is given in the last section IVTl 



II. EDM WITH ELECTRIC FIELD 



In our previous work 17| . we defined NEDM from the CP-odd electromagnetic form factor, 



F3, in the zero momentum transfer limit. In the actual calculation, however, it is not so easy 

7T 

to change the momentum transfer, since the momentum is quantized as p = n— on a finite 
spatial length of L. In the case of large p with n — 2, 3, • ■ ■ at small L, statistical errors 
are large, while a smaller p with n = 2, 3, • • • , which has a better signal, requires a larger 
lattice size L. In both cases, the calculation becomes more difficult for larger momentum at 
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n — 2, 3 • • • than for the smallest momentum at n = 1. In addition, the correct distribution 
of the topological charge is essentially important for the NEDM calculation. Since the 
width of the distribution of topological charge increases linearly with the volume, larger 
volume calculations require more statistics than the small volume ones, contrary to other 
observables. 

The difficulties for the extrapolation to the zero momentum transfer limit mentioned 
above are our motivations to consider a different method for the NEDM calculation with 
which we can avoid the momentum extrapolation. In this section we introduce our new 
approach for the lattice QCD calculation of NEDM. 



A. Formulation 

In ref. NEDM is defined through the energy change of the neutron state in the 
presence of an external electric field, similarly to the magnetic moment defined from that in 
the magnetic field. If a static and uniform electric field exists in a CP-violating system, the 
electric dipole moment (EDM) appears in the Hamiltonian as the interaction term between 
spin S of particle and electric field E: 

6H CP = d N (9)S-E + 0((E) 3 ) (7) 

where (1n{Q) represents the EDM. In order to extract EDM we consider the energy difference 
of nucleon states for opposite spins in the external electric field: 

mi(E) - m e _ g (E) = 2d N (6)S ■ E + 0((E) 3 ), (8) 

where m e ± g(E) denotes the energy of nucleon whose spin vector is ±S in the presence of 
the electric field E. Therefore we can extract cLn{0) from the nucleon propagators for 
two different spin stcitcs cit zero momentum only, avoiding difficult calculations at non-zero 
momenta. 

For small 6 we can expand d^(9) as 

d N {6) = d N 6 + 0(6 3 ). (9) 

We will check that higher order contributions at O(0 3 ) are negligibly small. Hereafter we 
represent djy as the leading order of EDM. 
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B. Methodology on the lattice 



A static and uniform electric field is represented by the spatial gauge potential as 

Ai(x) = Ef ud[d t (10) 

where Ef uchd is the constant electric field in Euclidean space. A non-zero NEDM could be 
detected from the oscillating behavior of the neutron propagator. Since NEDM is expected 
to be small, it is numerically very difficult to measure such a small oscillation. On the other 
hand, if we employ a static and uniform electric field in Minkovski space as 

Ai(x) = -iEf mkov t, (11) 

the oscillation turns into an exponential behavior, which is easier to measure. Therefore 
we introduce a static and uniform electric field in Minkovski space as an external field into 
lattice QCD, by replacing the spatial link variables as 

Ur\x) — ■* Ur\x-E^ m ) = e'^^Ur^x) , (12) 

where q e denotes the quark charge, 2/3 for up quark and —1/3 for down quark. Hereafter we 
suppress the superscript of the constant electric field Ei in Minkovski space for simplicity. 

An obvious problem here is that the Minkovski electric field E breaks the periodic bound- 
ary condition in the temporal direction: 

Ui(t + T,x) = Ui(t,x), (13) 
Ui(t + T, x ; Ei) = e^ E ' T Ui(t, x; Ei) ^ x; Ei) (14) 

where T is the size of the temporal direction. This generates an effective electric field, 
defined by E t (t) = Mt+V-Mt-i) ^ &g 



, Ej t — 2, 3, • ■ ■ ,T — 1 , 

m = { T2 ■ 

-Ei—y- t = 1, T 

Therefore the electric field is no more constant near the boundary between t — 1 and t = T. 
In order to avoid the effect of this non-uniform electric field to the EDM signal, we have to 
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take Ei as small as possible. In any case a small value of Ei is necessary to neglect 0((E) 3 ) 
terms in (JTj). 1 

In our calculation gauge configurations are generated by the usual lattice QCD action 
without E and 9. After inserting the electric field to gauge configurations we calculate quark 
propagators for flavor u and d separately, in addition to the normal one with Ei = 0, which 
is used to remove a fake signal at Ei = caused by statistical fluctuations. The total number 
of solvers for quark propagators is three for each configuration. From quark propagators we 
construct the nucleon propagator with the 9 term as 

(N a N p )g(E,t) = (N a (t)N a (0)e i6Q ) s (16) 

where (O) g represents the vacuum expectation value of O with E but without the 9 term. 
Here we use the re-weighting method with the complex weight factor e 18 ® . In order to obtain 
good signals, a large overlap of gauge ensembles between 9 = and 9 ^ as well as the 
correct distribution of the topological charge are required. Taking a small value of 9 as long 
as we get a signal helps for the large overlap, while we have to simply accumulate an enough 
number of configurations for the correct distribution of the topological charge. 

In the presence of the uniform and static electric field, the upper components of the 



nucleon propagator at zero spatial momentum take the following form for a, (3 = 1,2 
(N a N p )g{E,t) = Z' N \E 2 = E-E)[(l + A N {9,E 2 )a-E 



2l|: 



x exp ( - m e N (E 2 )t - M °' E >- B • Et)} + • • • , (17) 

2 J a/3 

where the EDM d^(9, E 2 ) and the spin-dependent amplitude A^(9, E 2 ) are odd in 9, while 
the spin-independent energy 2 m e (E 2 ) and an overall amplitude Z' N 6 (E 2 ) are even in 9. Here 
dots denote contributions from excited states. 

To extract EDM we construct the ratio of nucleon propagators between different spinor 
components. For E = (0, 0, E) we consider the following ratio: 

nnaive/p _ WVl) ((0, 0, S), *) _ 1 + A N (9, E 2 )E BFf > n <&F 9 F 3 )] flirt 

^ 3 {E ^ 6) - (N 2 N 2 ) e ((0AEU) ~ l-A N (9,E*)E eXp[ - dN6Et + 0{e E ^ )] ' (18) 



1 The electric field in Euclidean space smaller than Ei =2%/T also breaks the periodic boundary condition. 

2 The energy of the proton increases as t increase since the charged particle is accelerated in the uniform 
electric field. This effect is canceled in the ratio, which will be used to extract the signal of EDM. 
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where we use eq. ( IT7|) for the second equality. Similarly for E = (E,0,0) and (0, E, 0), we 
obtain 

(NxN 2 )e + (iV 2 #i) e + (iV 2 # 2 ) e 



R™ lve (E,t;- 



(iViiV 1 ) £ 
l + ^( 



(iViiV 2 ) e - (N2NJ0 + (iV 2 AT 2 ), 
£ 2 )£ 



-((£,0,0),t) 



(19) 



r> naive 

K 2 



(E,t;> 



T -^^eM-d N 0Et + O(e<E,0E% 
(NiNx)e + tiN^e - i<# 3 #i)* + (iV 2 iV 2 ) fl 



(iViiVi) ( 

1+Ajv( 



I - A 



- i{N l N 2 ) e + z(iV 2 iVi) e + (iV 2 iV 2 ) 



((0,£,0),t) (20) 



E 2 )E 



We can average over the ratio in three directions to increase statistics, if necessary. 

In order to remove the spurious contribution m|(0) — m 9 _g(0), which must vanish for 
infinite statistics, we consider a double ratio defined by 

~Rf aivc (E, t; 9 



In 



Ri(E, t; 9) 

Ri(E, t; 9) 
Ri(E,t + l;6 



#? aive (O,t;0) 
[miiBk) - m%0)\ - [mUEt) - m e _M] 
d N 9E + 0(9 3 E,9E 3 ). 



(21) 

(22) 
(23) 



We can improve the EDM signal further, removing the contribution at 9 = 0, which also 
vanish for infinite statistics, by a triple ratio as 



R 



(w/o 0=0) 



(E,t; 



Ri(E, t; 



1 + 9A l N (E 2 )E 



expldjyOEt], 



(24) 



Ri{E,t;9 = 0) 1 — 9A 1 N (E 2 )E 
where we used an expansion A^(9,E 2 ) = 9A l N (E 2 ) + 0(9 3 ), and we finally subtract the 
spurious contribution even in E by a quadruple ratio as 



R 



(w/o 0=0) 



(E,t; 



Rp /o9=0) (-E,t; 
R™ vc (E,t;9) R 



1VE (-E,t;9 



0) 



tii 



R™" c (E,t-9 = 0) 



M]v(E 2 )E 



exp[2div^t], 



(25) 



1-9A 1 N (E 2 )E, 

where the second equality tell us that this is indeed a triple ratio since E = contributions 
are canceled identically. We finally extract EDM from the exponential fit to R Eorr (E,t;9) 
over some time range, determined by the behavior of the effective EDM: 

Rf* T (E,t; 9) 



2d N 9E = In 



Rl orr (E,t 



1,2,3. 



(26) 



S 



III. SIMULATION DETAILS 



A. Simulation parameters 

In our study we employ gauge configurations generated by the RG improved gauge action 
at (3 = 2.6 in the quenched approximation, which corresponds to a" 1 = 1.902(50) GeV from 



the string tension a assuming a = (440 MeV )[22J. 



For the quark action, we employ the domain-wall fermion on a 16 3 x 32 lattice with the 
fifth length N s = 16 and the domain-wall height M = 1.8. These parameters are identical to 
those in the previous EDM form factor calculation. We however take a heavier quark mass, 
rrif = 0.12, which corresponds to m^a = 1.113(2) and mpsjniy ~ 0.88, than the one in the 
previous calculation, in order to reduce the computational cost, since our main motivation 
in this calculation is to see whether the EDM signal can be obtained by this method. As 
shown in the next section we have indeed obtained the EDM signal after accumulating 1000 
configurations at this heavier quark mass. 

We also investigate whether the EDM signal can be obtained by this method with the 
clover fermion. The EDM calculation with this fermion has the advantage that the compu- 
tational cost is roughly N s times smaller than the cost of the domain- wall fermion so that 
systematic studies such as volume or quark mass dependences can be performed more easily. 
Moreover we can employ the Nf = 2 and 2 + 1 flavor dynamical configurations already gener 



ated with the clover quark action at several sea quark masses and lattice spacings 23j, l24l. 125] 
in future studies. We calculate the EDM on the same 16 3 x 32 configurations, using the clover 
fermion with csw — 1-340, the tadpole improved value of the clover coefficient determined 
from 

csw=[J2 = t 1 - 0.8412/r 1 )" 1 / 4 . (27) 

X,fl<U 

In order to obtain a similar nucleon mass, we use the hopping parameter k = 0.1320, 
corresponding to m^a = 1.020(2) and mps/mv — 0.85. 

Since, as will be shown later, the EDM signal can be successfully obtained with the clover 
fermions, we investigate the volume dependence of the EDM signal using a 24 3 x 32 lattice. 
Furthermore the quark mass dependence of the EDM is calculated with this fermion on this 
larger volume. 

For the calculation of quark propagators we employ the smeared source of the form that 
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f(r) = Ae~ Br where r = \x — x src \ with the source point x src = (8, 8, 8) on 16 3 and (12, 12, 12) 
on 24 3 spatial lattice, after the Coulomb gauge fixing is applied to gauge configurations. We 
mainly take t src = 1 as the time slice of the smeared source. In order to check the effect of 
the non-uniform electric field near t — 1 and T, we also calculate the EDM with t src = 8, 
using; the clover fermion on a 24 3 x 32 lattice. Effective mass plots of nucleon in various 
cases are given in Fig. [TJ We observe the plateau at t > 7 for the domain-wall fermion and 
the clover fermion at heaviest quark mass, while plateau appears at t > 6 for the clover 
fermion at lighter quark masses. 

In our calculation we mainly take E = (0, 0, E) with E = ±0.004. As exceptions, 
E = ±0.002 is employed on a 16 3 x 32 lattice with the domain-wall fermion to investigate 
the E dependence of the EDM signal, and (E, 0, 0) and (0, E, 0) are used on a 24 3 x 32 lattice 
with the clover fermion at heaviest quark mass to check the consistency and to increase 
statistics. Although we can easily change the value of 9 by reweighting, we fix 9 = 0.1 in 
our calculation, except 9 = 0.05 and 0.2 on a 16 3 x 32 lattice with the domain-wall fermion 
to investigate the 9 dependence of the EDM signal. 

Parameters of fermion actions in various cases are summarized in Table HI 



B. Topological charge 



The topological charge using the 0(a 2 ) improved definition 26J is measured on each con- 
figuration after 20 cooling steps. 

On a 16 3 x 32 lattice we accumulate 1000 configurations. In Fig. [2] we present the 
histogram of the topological charge, which is consistent with gaussian distribution. The 
symmetry of the distribution is measured by the average of Q, which is consistent with zero 
within error: (Q) = —0.002(97). If the gaussian distribution is assumed, its width a is 
given by (Q 2 ) = 9.37(44). On this lattice size 1000 configurations seem enough to give a 
reasonable distribution of the topological charge. 

On a larger volume of a 24 3 x 32 lattice, we accumulate nearly 2000 configurations since 
(Q 2 ), thus the width of the distribution of Q, increases linearly in volume. In Fig. [3], we 
show the histogram of Q, which looks reasonable, namely sufficiently symmetric and close 
to gaussian. We find (Q) = 0.15(13) and (Q 2 ) = 33.6(1.1). 
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IV. EDM SIGNAL AND SYSTEM ATICS 



In this section, we show numerical results for nucleon EDM signals with the external elec- 
tric field method. We investigate several systematics of the EDM signal such as dependences 
on the fermion action, the volume, E, 9, t src and the direction of E. 

A. Comparison between domain-wall and clover fermions 

We first consider the case of the domain-wall fermion on a 16 3 x 24 lattice. In Fig. H] 
we plot the double ratio R^(E,t;9) clS db function of t at (E,6) = (±4.0 x 10" 3 ,0.1) and 
(E,6) = (±4.0 x 10~ 3 ,# = 0), for both neutron and proton. The star symbols in Fig. HI 
representing the time dependence of R^(±E, t; 9 — 0), are consistent with unity within errors 
at both ±E. This confirms the expected behavior that the exponential part of R^{E,t]9) 
vanishes at 9 = 0. For non-zero 9, on the other hand, deviations of Rs(E,t;9) from unity 
show up beyond errors and they increases as t increases. Moreover the sign of deviations 
depends on the sign of E. All these behaviors of R3 are consistent with the fact that non- 
zero value of EDM exists. In Fig. [5] we plot time dependence of R^^ 09-0 ^^, t; 9), defined 
in eq. (124")) . for which contributions at 9 = due to finite statistics are removed. The E 
dependence of signals become more visible after the removal of 9 = contributions. In 
addition it is noted that the EDM signal of proton has an opposite sign to that of neutron. 

Applying the same analysis as above to the case of the clover fermion on a 16 3 x 32 lattice, 
we obtain a similar behavior for R3 and j^/ o9 -°\ Therefore we do not present them here. 
Instead the effective mass of i?3° rr , defined in eq. (125]) . is plotted as a function of t in Fig. [6l 
for both domain-wall and clover fermions. It is interesting to see that the time dependences 
of the effective mass for the two fermions are very similar. Moreover, for both fermions, 
we observe plateau around 6 < t < 12, whose values are non-zero beyond errors. Clearly 
the EDM signal for proton has an opposite sign to that for neutron, as suggested by the 
behavior of ji^/ oe -°\ 

Let us conclude this subsection. Using the external electric field method, we obtain 
the EDM signal for both neutron and proton, with both domain-wall and clover fermions. 
This suggests that the chiral property of the fermion action does not play a crucial role 
to obtain the EDM signal with this method. Note however that the quark mass employed 
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in this investigation is rather heavy. Therefore there is a possibility that some qualitative 
difference between two fermion formulations may show up at lighter quark mass where the 
chiral symmetry becomes important. In the remaining of this paper, we mainly employ the 
clover fermion formulation. 

B. Volume dependence 

We investigate the volume dependence of the EDM signal on a 24 3 x 32 lattice with the 
clover fermion at the heaviest quark mass. Here the physical spatial volume is increased to 
2.4 3 fm 3 from 1.6 3 fm 3 . Our main concern is whether the nonzero value of the EDM signal 
obtained in the previous subsection persists as the volume increases. 

In Fig. Owe compare the effective mass plot of R^ orr (E, t; 9) at 9 = 0.1, E = 0.004 in the 
larger volume with that in the smaller volume. It is clear that the EDM signal remains non- 
zero in the larger volume. Results in both volumes are consistent with each other within 
large errors. We can conclude that the EDM signal obtained with this method does not 
vanish in both volumes. 

C. Boundary effect of the electric field 

The electric field in our method breaks periodicity in the time direction, leading to a large 
non- uniformity near the boundary between t = 1 and t = T. Since we put a source at t = 1, 
the EDM signal may be affected by the non-uniform electric field. In order to investigate 
how the EDM signal is affected by this boundary effect, we repeat the EDM calculation on 
a 24 3 x 32 with the clover fermion at the heaviest quark mass, moving the source point to 
the different time slice but keeping other conditions fixed. 

In the previous calculation at t src = 1, we observed that the plateau seems to exist at 
t > 8. Since this indicates that the effect of boundary may be small at t — 8, we take a 
new source point at t src = 8. If we need a minimum plateau length of 5 for a reliable fit, we 
wonder be using a plateau at t — 15 — 19 for t src = 8. Since the time slice t = 19 or 20 is 
largely separated from the boundary at t = T = 32, the boundary effect to the plateau as a 
whole is expected to be small. Therefore t scr = 8 is a reasonable choice. 

In Fig. [8] we compare the time dependence of R^^° e ~°\±E,9,t) for two different source 
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points, t src = 1 and t src = 8. We clearly observe a different time dependence of ^( w / o6, -°) f or 
two sources at small time slices, t — t src < 4. We think that large deviations of ^( w / 061 - ) f r0 m 
unity at t — t src < 4 for the t src = 1 case is an effect of the large non-uniform electric field 
near the boundary between t — 1 and t = T. On the other hand, the deviation of j^/ o9 -°} 
from unity becomes visible around t — t src ~ 4 for the case of t src = 8. Since the plateau of 
the nucleon effective mass appears around t — t STC ~ 5 — 6, contributions from excited states 
to /£( w / oe -°) become small and the nucleon state dominates around this range of t in the 
case of t src = 8. In Fig. [9]we plot the effective mass of Rf TT (E, 9 = 0, t) for the t src = 8 case, 
together with that for the t src = 1 case. We notice that the plateau starts around t — t src = 5 
for the t src = 8 case. For the t src = 1 case, on the other hand, the values of effective mass of 
R^ orr (E, 6 — 0, t) around t — t src = 4 — 6 seems smaller than the plateau of the t src = 8 case, 
suggesting that the boundary effects, observed in ^( w / oe -°) a ^ sma ll t — t src , still remain in 
the effective mass around t — t STC = 4 — 6. Therefore, to avoid possible contaminations from 
the boundary effect, we take sufficiently large separations such that t — t STC = 8 — 11 for the 
fit of R c 3 orr (E, 9 = 0, t) in the case of t src = 1. 

An important lesson here is that we should take the starting point of the fitting range as 
far from the source as possible, if the source is placed near the boundary such as t src = 1. 
This caution should be applied to all other data obtained with t STC = 1 . 

Fitting with i?3° rr exponentially in 5 < t — t STC < 9 with t SIC = 8, we obtain 

\ -0.025(8) e ■ fm (Neutron) 
d N = < , (28) 

0.024(11) e ■ fm (Proton) 

while for the t STC = 1 case we have 

f -0.030(8) e • fm (Neutron) 
d N = < (29) 

0.036(11) e-fm (Proton) 

with t — t src G [7, 11] as the fitting range. Two results are consistent with each other within 
large statistical errors. Similarly, on a 16 3 x 32 lattice, we obtain 

-0.021(11) e-fm (Neutron) 
d N = < (30) 
0.026(13) e-fm (Proton) 



for the clover fermion and 



-0.017(8) e • fm (Neutron) 
d N = < (31) 
0.020(10) e-fm (Proton) 
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for the domain- wall fermion. The fitting range is t — t src G [6, 11] with t src = 1 for both 
fermions. These values, summarized in Table [Til have the same sign and a similar order of 
magnitude to the EDM form factor previously obtained on a 16 3 x 32 lattice with the domain- 
wall fermion with the form factor method, which is given by F 3 (q 2 ~ 0.58GeV 2 )/mjv = 
—0.024(5) e-fm for neutron and 0.021(6) e-fm for proton 17|. These agreements of sign and 
magnitude between the two methods support that the viability of this method explored in 
this paper. 



D. E and 9 dependence 

In Fig.nUlwe plot values of EDM as a function of E for neutron(upper) and proton (lower) 
at 9 = 0.1. Observing the expected linear dependence on E for both cases, we conclude that 
0(E 3 ) contributions in (jSJ) are negligible. Fig. [TT] shows djsi{9) in lattice unit as a function 
of 6 at E — 0.004, assuming the linear E dependence of the fitted EDM signal. We again 
confirm that the linearity in 9 is good and thus 0(9 3 ) contributions in ([9]) are reasonably 
small. 

We concluded that our choices of {E,9) = (0.004,0.1) are small enough to ensure linear 
dependences of the EDM signal on both E and 9, which we assume in the analysis in the 
rest of this paper. 



E. Average over the electric field 

Averaging over three directions of the electric field is not so useful in quenched simulation. 
This way of increasing statistics, however, may become important in full QCD case since 
the number of full QCD configurations is limited. In this subsection we investigate the 
effectiveness of this method and the related question of the independence of the EDM signal 
on the direction of the electric field. 

Using eq. (Ell), eq. (EEH]) and eq. !j20l) for E = (0,0, E), (£,0,0), and (0,£,0), we obtain 
Ri as a function of E for each £ on a 24 3 x 32 lattice with the clover fermion at heaviest 
quark mass. In Fig. [TJl Ri shows similar time dependences for all % . EDM signals, given 
in Fig. [131 are also comparable in the similar time range among different directions. We 
confirm the consistency among extraction of the EDM signal from three different directions 
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using the formulae in eqs. (fl8]) - (l20l) . 

We now consider the average over 3 directions. In Fig. [TUthe effective mass of the average, 
R COTT (E,9,t) = J2i Rf TT (E, 6>, t) is plotted as a function of t. Fitting it exponentially at 
7 < t — t src < 11, we obtain 

f -0.0276(72) e • fm (Neutron) 
d N = < . (32) 

[ 0.0278(87) e • fm (Proton) 

Although errors are reduced in the effective mass, the reduction in is much smaller than 
1 / \f?>. We conclude that the error reduction by this averaging is limited, due to the possible 
correlation among R i= i i2y 3(E,9,t), 

V. QUARK MASS DEPENDENCE 

In this section we study the quark mass dependence of EDM using the clover fermion on 
a 24 3 x 32 lattice. 



A. Quenched effects 

It is well known in full QCD that EDM generated by the 9 term must vanish in the chiral 
limit. This can be seen from the fact that the CP- violation Lagrangian after an appropriate 
chiral rotation Q], 

5C C p = iOrh ^ m = ( 5^ m i 1 ) ' ( 33 ) 

i=u,d,s i=u,d,s 



vanishes in the massless limit of any quarks. (See 20J for more detailed argument on this 
property.) 

In quenched QCD, however, this argument fails since the 9 parameter can not be trans- 
lated to the above form in the absence of the chiral anomaly, which requires the quark 
determinant. Therefore CP- violating observables generated by the 9 term may remain non- 



zero in the zero quark mass limit. Indeed, as discussed in [20], zero modes of the quark Dirac 
operator can generate CP-odd contributions even in the massless limit. It is not so easy, 
however, to determine the explicit quark mass dependence of the EDM from the general 
argument in quenched QCD. 
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Recently, from the numerical simulation of the instanton liquid model 27[, the 
dependence for NEDM has been reported near the chiral limit of quenched QCD. The 
partially quenched chiral perturbation theory 28|, on the other hand, has suggested the 
1/m 3 behavior in the finite volume of L 3 at fixed sea quark mass m sea such that 



,P.Q.ChPT 

a N ' sTs" 

from the leading contribution of one-loop graphs. 



fn i 



(34) 



B. Quark mass dependence of EDM 

We calculate EDM at three different quark masses with the clover fermion on a 24 3 x 32 
lattice. In Fig. [TS] and Fig. [T6]we plot the effective mass of R corr (E, 9, t) = £ 3 =1 #9° rr (£, 9, t) 
as a function of t at two lighter quark mass with t STC = 1. Signals become a little noisier and 
less stable as the quark mass decreases. Fitting data at t — t src G [7, 10] for the three quark 
masses, we obtain the quark mass dependence of EDM for neutron and proton as shown in 
Fig. [13 and Table IIHI Compared with the current algebra result, —0.0036 e-fm {?], Q] also 
shown in the top of Fig. [T71 our quenched NEDM are about 10 times larger. Moreover our 
results suggest that EDM does not vanish in the chiral limit for both neutron and proton. 
We consider that the larger value of NEDM we focus is partly due to the quenched effect. 
Because of large statistical errors, we can not distinguish the functional form of the mass 
dependence of EDM, whether it stays constant or diverges in the chiral limit. 



C. Quark mass dependence of the CP-odd phase factor 

In addition to the EDM, using the clover fermion, we calculate a simpler quantity flf, 



the CP-odd phase factor of the nucleon propagator, defined in Ref. 



Uj as 



(N(p } t)N(p } 0)Q) = l^v|V^4^5. (35) 

Since the CP-odd phase factor arises from CP-violation effects of the 9 term, would 
vanish in the chiral limit of full QCD. In quenched QCD, however, this quantity also may 
remain non-zero in chiral limit because of the same reason as the EDM. 

In Fig. [THl we show the time dependence of the nucleon propagator at the next leading in 
9, — tr[(N(0, t)N(0, 0)<5)-f ] (left), and effective masses of the leading nucleon propagator in 
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9, tr[(A r (0, t)N(0, 0))^-^}, as well as the next leading one (right) at 3 quark masses. Since 
effective mass plots show the agreement of masses between two propagators around t = 10, 
we extract fa by fitting tr[(iV(0, t)N(6, 0)Q)^] at 9 < t < 12 in the form of ([35}, where 
\Z^\ 2 and m^v have been fixed from the leading propagator. 

The quark mass dependence of fa is given in Fig. [H and Table DTD It is noted that 
errors of fa are much smaller than those of EDM. The top of Fig. [19] shows that fa does 
not vanish in chiral limit and moreover it seems to diverge as l/m q in this limit. To see this 
behavior more clearly, we plot /jy multiplied by the quark mass m q = (k _1 — n~ l )/2 as a 
function of m q in the bottom of Fig. [19j The fact that /^m 9 seems almost constant at this 
range of the quark mass suggests that fa may diverge as l/m q in the chiral limit. It may 
be interesting to confirm this behavior of by some theoretical considerations. 

VI. SUMMARY AND DISCUSSION 

In this paper, we have investigated the viability of an old idea for calculating the nucleon 
EDM by introducing a uniform and static electric field. In this setup the nucleon EDM ap- 
pears directly in the energy difference between spin-up and spin-down states of the nucleon. 
To introduce the complex 9 term into lattice QCD calculations, we used the reweighting 
technique with the factor e ie ® . We have demonstrated that this reweighting method indeed 
works as long as 9 is small enough, by calculating the nucleon EDM in quenched QCD on 
a 16 3 x 32 lattice at a relatively heavy quark mass. We found that the quality of signals is 
not very sensitive to lattice fermion formulations employed, domain-wall fermion and clover 
fermion in our study. Using the clover fermion on a 24 3 x 32 lattice, we investigated the 
effect of non-uniformity of our electric field induced at the boundary in time direction. Even 
if the source point of nucleon is placed near the boundary, the effect to the nucleon EDM 
disappears for large enough t, while the effect becomes smaller even at small t if the source 
is placed away from the boundary. We also found that the finite size effect to EDM is not 
so large: results between (1.6 fm) 3 and (2.4 fm) 3 boxes agree within errors. 

We investigated the quark mass dependence of the nucleon EDM and the CP-odd phase 
factor in quenched approximation on a larger volume with the clover fermion. Both 
quantities do not seem to vanish in the chiral limit, in contrast to full QCD where effects 
of the 9 term disappear for a massless quark. Therefore non-vanishing behaviors of EDM 
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and flf are purely quenching effects. In particular, seems to diverge as 0(l/m q ) in 
chiral limit. It is, however, difficult to determine precise quark mass dependences of these 
quantities in quenched QCD, due to larger statistical errors. 

This work shows that the external electric field method is simple and straightforward for 
the determination of the EDM in lattice QCD. In particular, the success with clover fermion 
in this method is significant for applications to full QCD simulations. We are currently 
carrying out the EDM calculation using Nf = 2 dynamical clover configurations generated 
by the CP-PACS collaboration [19 1. 
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APPENDIX A: ELECTRIC POL ARIZ ABILITY OF NEUTRON 



In this appendix we discuss the electric polarizability of the neutron. This observable 
can also be obtained by the external field method employed in our calculation, as has been 



done in refs. 



29 



301 ] . We compare our results with theirs. 



1. Definition 

The electric polarizability on is defined as the coefficient of the E 2 term in the expansion 
of the E dependent nucleon mass rriN^E): 

Am N {E) = m N (E) - m N (0) = -^Tra^Ke -1 ^ 2 ^) 2 . (Al) 

which is measured by Compton scattering experiments. Note that the electric field E here 
is dimensionless. A recent Compton scattering experiment gives 

o^ p = (1.16 ±0.15) x 10~ 3 fm 3 (A2) 
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for the neutron 



311 ] . In the lattice calculation the effective mass shift is calculated by 

(NN) (E, t) 



r N (E,t) 
Am N (E) = In 



(NN){0,t)' 
r N (E,t) 



(A3) 
(A4) 



r N (E,t + l) 

where (NN)(E,t) denotes the nucleon propagator in the presence of the constant electric 
field E without reweighting e 10 ®. In order to remove spurious contributions odd in E from 
the effective mass shift, we take an average over E and —E, by replacing 

1 



r N (E) 



■(r N (E,t) + r N (-E,t)) 



(A5) 



in eq-fTAH). 



2. Numerical results on a 16 3 x 32 lattice 



Our lattice setup for the calculation of the electric polarizability is same as the one 
employed for the NEDM calculation in sec. IIV Al In particular, the real electric field E = 
(0,0, E) in Minkovski space is introduced by the replacement of eq. ffT2j) . Although the 
periodicity in time direction is broken by this electric field, the boundary conditions for the 
fermion are periodic in both time and spatial directions on a 16 3 x 32 lattice. We employ the 
domain-wall fermion at E = 4 x 1CT 3 and E = 2 x 1CT 3 . As a comparison we also employ 
the clover fermion at E = 4 x 10 -3 . 

In the top of Fig. [20] we show the effective mass plot of in eq. (1A4|) for domain-wall 
and clover fermions on same configurations. We observe the plateau starting around t = 7 
for the clover fermion and around t = 10 for the domain-wall fermion. From the exponential 
fit of rjv(t) at 9 < t < 14, we obtain Am^, whose values are given in Table HVl 

In the bottom of Fig. [20] we present the E dependence of the mass shift 2ArriN for 
the domain-wall fermion. By fitting data with — 4:Tra^(e 2 a i )~ 1 E 2 , we obtain the electric 
polarizability for neutron : 

a N = 1.32(2) x 10" 4 fm 3 (A6) 

in the unit of e 2 (47r)~ 1 a 3 ~ 0.73 x 10~ 5 fm 3 with the fine-structure constant a = e 2 /(A7r) = 
1/137. 

This value, obtained in quenched QCD at a ~ 0.1 fm and mps/my ~ 0.88. is 1/10 times 
smaller than the experimental value a e ^ p = 1.16(15) x 10~ 3 fm 3 , but the sign of agrees. 
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3. Results on 24 3 x 32 with two different source points 

We also calculate the electric polarizability of the neutron on a larger volume, 24 3 x 32, 
using the clover fermion at k = 0.1320. As in sec. IIVCI we employ two different source 
points, t src = 1 and t src = 8, to investigate the effect of the gap in E at the boundary to the 
electric polarizability. 

In Fig. |2T1 we present the effective mass shift, Am^a, for both t STC = 1 and t src = 8. 
Compared with results on the smaller volume in sec. IA 21 plateaus seem to appear at very 
large t for both sources or even Arriba may not reach the plateau at t < 16. Even though an 
identification of plateaus is less reliable on the larger volume, we fit data exponentially in t 
at 13 < t < 16 and give values of Arriba in Table HVl As seen in the table, the magnitude of 
fitted values is larger than the value on the smaller volume. We think that this discrepancy 
is mainly caused by contaminations from excited states on the larger volume. We need larger 
time separations to extract the ground state contribution unambiguously. We also observe 
large differences in the effective mass at small t between t src = 1 and t src = 8. This indicates 
that the electric polarizability is quite sensitive to the boundary effect. 

In Fig. [22] we plot the effective mass shift at different quark mass after taken average over 
three directions of electric field with t src = 1 on 24 3 x 32. We observe that the time behavior 
is not so different with each other, and therefore its value will not depend on the quark mass 
strongly. Fig. [23] and Table [V] shows the converted results to electric polarizability using 
fitting data of Amjv in each k. In these heavier masses, the results seems to be constant 
for square of pion mass, though statistic errors are still large. Therefore more statistics are 
probably needed to give a precise value of the neutron electric polarizability in the chiral 
limit. 



4. Comparison with previous calculations 



29 



30|: 



As a test of our method, we use same lattice parameters as in previous calculations 
Accumulating 40 quenched configurations generated by the plaquette action at f3 — 6.0 ( 
a ~ 0.1 fm ) on a 24 4 lattice, we calculate the electric polarizability by the Wilson fermion 
action at n = 0.1515, which is the heaviest quark mass in 30]]. With the periodic boundary 
condition in spatial directions but the Dirichlet boundary condition in the time direction, 
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the nucleon propagator is calculated for a point source at t = 1 and a point sink at t. 
The electric field is introduced into all spatial link variables in the expanded form: 



U 3 {x) — > e i(lEt U z {x) ~ (1 + iqEt)U 3 {x), (A7) 

where we use an electric field in Euclidean space, which corresponds to the imaginary value 
in Minkovski space. Therefore the E dependence of the mass shift Amjv is given by 

1 1 

Am N (iE) = --(A-Ka N ){ie- l a- 2 E) 2 = -(A-Ka N )e- 2 a- A E 2 (A8) 



with the electric polarizability a^. As in |30|], we employ E = ±1.08 x 1CT 3 , ±4.32 x 
10 -3 , ±8.64x 10~ 3 in the actual calculation. Note that the periodicity of spatial link variables 
in the time direction is explicitly violated partly due to the fact that E ^ 2ir/L and partly 
due to the expansion (1A7I) . 

Fig. [241 shows the effective mass shift for neutron in eq. (1A4|) at \E\ = 1.08 x 10 -3 . Our 
data in Fig. [24] roughly agree with filled circle symbols in Fig. 6 of 30j . Unfortunately a 
candidate for a possible plateau appears only atl5<£<19. Assuming that this is indeed 
a real plateau, we fit Am at exponentially in t at 15 < t < 19 and gives values at each E in 
Table El 

In Fig. [25] we plot the E dependence of mass shift Am^. By fitting data with 
|(47rajv)e _2 a _3 .E 2 , we obtain a coefficient a at, the value of electric polarizability: 

Kbrichict = -8-5(8) x 10" 4 fm 3 . (A9) 



This value agree with the value in [30j, [aArjomchiet = —7.9(5) x 10~ 4 fm 3 , within about 
one-sigma error. Surprisingly the sign of this result is opposite to the result flA6j) obtained 
by the real electric field in Minkovski space and to the experimental value in eq. (1A2[) 3 . In 
addition we confirm that the negative value of is obtained even if we use the real electric 
field in Minkovski space in the Dirichlet boundary condition. Therefore the wrong sign of oin 
in this case is not caused by the way of introducing the electric field (Euclid or Minkovski) 



In 
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301 ] it has been claimed that the electric field inserted as in eq. (|Alj) is real so that their results of 
the electric polarizability have the same sign as the experimental value in eq. (|A2|) . However, as shown 
here, the electric field introduced by eq. (|A1[) is real in Euclid space and it becomes pure imaginary in 



Minkovski space. Therefore electric polarizabilities in 29|, [30( are opposite in sign to the experimental 
value. 
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but is related to the boundary condition in the time direction. We think that T = 24 is 
too short to suppress contributions from excited states to ajv- In order to obtain a reliable 
estimate for a^, one should investigate dependences of results on the lattice set-up such as 
the boundary conditions, the source point or the way of introducing the electric field. We 
leave these studies in future investigations. 
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FIG. 1: The effective mass plot for nucleon with domain-wall fermion at m-f = 0.12 (top-left), 
clover fermion at k = 0.1320 (top-right) on a 16 3 x 32 lattice and clover fermion at various quark 
masses (bottom) on a 24 3 x 32. 
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FIG. 2: (top) Histogram of topological charge improved by 0(a?) after 20 cooling steps. The solid 
line denotes the expected gaussian distribution from a = 

V(Q 2 ) ~ (Q) 2 - (bottom) The topological 

charge in each configuration. 
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FIG. 3: (top) Histogram of the topological charge in 24 3 x 32 lattice and (bottom) the topological 
charge in each configuration as shown in Fig. [2j 
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FIG. 4: The time behavior of R 3 (E,t;0) in E = ±0.004, = 0.1 and E = 0.004, 9 = with 
domain-wall fermion. (Top) neutron case, (bottom) proton case. 
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FIG. 5: The time behavior of R 



( ^ /o8 ~ 0) (E,t;e) inE = ±0.004, 9 = 0.1 with domain-wall fermion. 



(Top) neutron case, (bottom) proton case. 
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FIG. 6: The effective mass plot of Rf rr (E, t; 9) in E = 0.004, 9 = 0.1 with domain-wall and clover 
fermions in 16 3 x 32 lattice, (top) neutron case, (bottom) proton case. 
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FIG. 7: The same figure as Fig. [6] for clover fermion in 16 3 x 32 and 24 3 x 32 lattice, (top) neutron 
case, (bottom) proton case. 
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FIG. 9: These figures show the effective mass plot of the spin ratio, Rf vv (E, t; 8), in both t src = 1 
and t src = 8 at 9 = 0.1 and k = 0.1320. 



32 




0.001 0.002 0.003 0.004 

E 



FIG. 10: The E dependence for the fitting results of Rf"(E, t; 6) in the range of [7, 12] at 6 = 0.1 
with domain-wall fermion (top) for neutron, (bottom) for proton. The solid line denotes a result 
of linear fit. 
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E = 0.004 with domain-wall fermion (top) for neutron, (bottom) for proton. The solid line denotes 
a result of linear fit. 
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FIG. 12: The time dependence of Ri(E,6,t) for three direction of electric field, which 
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electric field, (top) neutron case, (bottom) proton case. 
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case, (bottom) proton case. 
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FIG. 14: These figures show the effective mass plot of the average of Rf orr (E,9,t) for three direc- 
tions of electric field with k = 0.1320. (top) neutron case, (bottom) proton case. 
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FIG. 15: The same figure as Fig. with k = 0.1330. (top) neutron case, (bottom) proton case. 
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FIG. 16: The same figure as Fig. [Ill with k = 0.1340. (top) neutron case, (bottom) proton case. 
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FIG. 17: The mass dependence of EDM factor with clover fermion. In top figure the star symbol 
shows the prediction from current algebra in (top) neutron case, (bottom) proton case. 
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FIG. 18: In the left figure we show that the time dependence of nucleon propagator in next leading 
order of at each quark mass: —ti(NNQ)^-. The straight line represents fitting function with 
f(x) = Ae~ mNt . In right figure we show that the comparison with exponent of nucleon propagator 
between the leading(open circles) and the next to leading order (solid circles). 
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FIG. 19: The mass dependence of CP-odd phase factor with clover fermion. The top figure presents 
/jy at each quark mass, and the bottom figure presents the /jy multiplied by quark mass. 
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FIG. 20: (Top) The effective mass shift plot as Fig.[25]in periodic boundary condition with domain- 
wall and clover fermion at E = 4.0 x 10~ 3 . (Bottom) The E dependence of the mass shift of the 
above results from exponentially fitting in the range of [9, 14] with domain- wall fermion. The 
solid line presents fit results with the function f(E) = aE 2 . 
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FIG. 21: This figure shows the comparison with different source point t STC . We plot the effective 
mass shift as Fig. [25] in periodic boundary condition with clover fermion at k = 0.1320 and 
E = 4.0 x 10~ 3 in large lattice size 24 3 x 32. 
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FIG. 22: We plot the effective mass shift as Fig. [21] with several quark masses. 
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FIG. 23: This figure shows the mass dependence of electric polarizability of neutron with same 
parameter as Fig. | 
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FIG. 24: The effective mass shift plot for neutron with Wilson fermion at E = 1.08 x 10 3 and 
k = 0.1515, The initial time of our data is set up to be same as [3C|. 
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TABLES 



TABLE I: Table for lattice parameters. The column of (A, B) denotes the smearing source param- 



eter in the exponent 


Fermion j3 L 3 x T x N s 


M a- 1 [GeV] 


m q 


(A, B) mps/my m^a 


Domain-wall 2.6 16 3 x 32 x 16 


1.8 1.902(50) 


0.12 


(1.28,0.40) 0.8781(4) 1.1130(15) 


Fermion (3 L 3 x T 


c sw a" 1 [GeV] 


K 


{A, B) mps/mv m^a 


Clover 2.6 16 3 x 32 


1.340 1.902(50) 


0.1320 


(1.55,0.24) 0.8508(5) 1.0202(17) 


Clover 2.6 24 3 x 32 


1.340 1.902(50) 


0.1320 


(1.55,0.35) 0.8494(1) 1.0186(9) 






0.1330 


(1.55,0.31) 0.8026(2) 0.9058(14) 






0.1340 


(1.55,0.27) 0.7253(2) 0.7843(16) 






= 0.1359(1) 
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TABLE II: Table for EDM results in some lattice parameters 



fermion m^a lattice size source point fitting range cZjv (Neutron) (Proton) 

domain-wall 1.1130(15) 16 3 x 32 t SIC = 1 t- t STC £ [6, 11] -0.0170(79) 0.0196(95) 

clover 1.0202(17) 16 3 x 32 t STC = 1 t — t STC E [6, 11] —0.0205(104) 0.0256(125) 

clover 1.0186(9) 24 3 x 32 t STC = 1 t — t src G [7, 11] -0.0304(78) 0.0361(111) 

clover 1.0200(9) 24 3 x 32 t STC = 8 t - t SIC e [5,9] -0.0246(83) 0.0237(112) 
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TABLE III: The mass dependence of EDM factor from exponential fit in the range 8 < t < 12 
for R cori (E , , t) which is the average over indices z = 1,2,3 and CP-odd phase factor in the next 
leading term of nucleon propagator. 



Neutron 




Proton 




fh 


fN m q 


K fit 


(In [e • fm] 


fit 


(In [e • fm] 







0.1320 -0.000212(48) -0.0276(72) 0.000214(67) 0.0278(87) -0.1075(80) -0.0117(8) 
0.1330 -0.000276(67) -0.0359(87) 0.000271(95) 0.0353(123) -0.1653(111) -0.0133(9) 
0.1340 -0.000300(97) -0.0391(125) 0.000300(143) 0.0390(187) -0.2738(152) -0.0143(8) 
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TABLE IV: Summary of the fitting results of mass shift of neutron with different boundary condi- 



tion and fermions. 



gauge action 


mass 


lattice size 


B.C. 


^src 




E 


Arriba 


Domain-wall fermion 


RG Iwasaki (3 = 2.6 


m f = 0.12 


16 3 x 32 


Periodic 


^src 


1 


Real, 0.002 


-0.0000375(44) 














Real, 0.004 


-0.000157(18) 


Clover fermion 


RG Iwasaki (3 = 2.6 


k = 0.1320 


16 3 x 32 


Periodic 


^src — 


1 


Real, 0.004 


-0.000155(20) 






24 3 x 32 


Periodic 


^src — 


1 


Real, 0.004 


-0.000265(22) 










^src — 


8 


Real, 0.004 


-0.000356(50) 


Wilson fermion 


Plaquette (3 = 6.0 


k = 0.1515 


24 3 x 24 


Dirichlet t src = 


1 Imag, 0.00108 


-0.000069(2) 














Imag, 0.00432 


-0.00107(18) 














Imag, 0.00864 


-0.00435(65) 
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TABLE V: Summary of the fitting results of polarizability of neutron with clover fermion action 
at three different quark masses after average over three directions of electric field. 

gauge action lattice size B.C. mass i src E a7\r(fm~ 3 ) 

RG Iwasaki (3 = 2.6 24? x 32 Periodic k = 0.1320 t src = 1 Real, 0.004 0.000227(14) 

k = 0.1330 t STC = 1 Real, 0.004 0.000226(26) 
K = 0.1340 t STC = 1 Real, 0.004 0.000228(63) 
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